Exact Solution for Vortex Dynamics in Temperature Quenches of Two-Dimensional Superfluids 



Andrew Forrester, Han-Ching Chu, and Gary A. Williams 

Department of Physics and Astronomy, University of California, Los Angeles, CA 90095 

(Dated: November 12, 2012) 

An exact analytic solution for the dynamics of vortex pairs is obtained for rapid temperature quenches of 
a superfluid film starting from the line of critical points below the critical temperature Tkt ■ An approximate 
solution for quenches at and above above Tkt provides insights into the origin of logarithmic transients in the 
vortex decay, and is in general agreement with recent simulations of the quenched XY model. These results 
confirm that there is no "creation" of vortices whose density increases with the quench rate as predicted by the 
Kibble-Zurek theory, but only monotonic decay of the thermal vortices already present at the initial temperature. 
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Although the phase-ordering kinetics of temperature- 
quenched thermodynamic systems have been studied for 
many decades QJJ] , there are only a few exact results for the 
dynamics of the recovery to equilibrium Si . For many sys- 
tems, progress in the field has been made by asserting that 
dynamic scaling should apply to a quenched system, that for 
a system with non-conserved order parameter the dynamics 
will be characterized by a growing length scale £(t) = ^t 1 /*, 
where z « 2 is the dynamical exponent of model A in the 
classification of Halperin and Hohenberg [4] and t is the time 
from the quench to low temperature. Scaling holds if solu- 
tions involving a length scale r only depend on the ratio r/£. 
The growing length scale characterizes the domain growth of 
the topological defects of the order parameter as the system 
becomes completely ordered at long times. A phenomenolog- 
ical dimensional argument |li|5t] is commonly used to predict 
the time dependence of the decaying defect density: 



P (t) oc f n oc r n ' z 



(1) 



where n is the number of components of the order parameter. 
Superfluids are in the n = 2 universality class, where the de- 
fects are quantized vortices, so if z — 2 then dynamic scaling 
predicts a 1/t decay of the vortex density. Computer simula- 
tions of spin systems with varying n gave general agreement 
with Eq. (Q~|), though only at long time and length scales [6]. 

However, a problematic case for dynamic scaling has been 
two-dimensional superfluids, where the defects are the vortex 
pairs of the Kosterlitz-Thouless theory |7] characterizing the 
equilibrium phase transition occurring at the critical temper- 
ature Tkt- The above arguments would give a vortex den- 
sity decaying as i -1 for quenches from well above Tkt to 
very low temperatures. Simulations of the XY model 181 41011 . 
however, showed the initial vortex decay to be considerably 
slower than this, and then only at very long times finally ap- 
proached the predicted exponent of —1. The behavior could 
be modeled as a \nt/t variation, but this requires altering the 
dynamic length scale to vary as (t/hit) 1 ^ for initial tem- 
peratures above Tkt- This change has been cited as a 
breakdown of dynamic scaling, though others lfl2ll find such 
a sudden change in the dynamic scale still fully consistent 
with scaling. To further complicate the issue, numerical so- 
lutions of the Fokker-Planck equation for vortex pairs carried 



out by two of us [13] for quenches starting from Tkt showed 
a rather different logarithmic form for the decay, as 1/ (t In t), 
and quenches from initial temperatures below Tkt showed a 
temperature-dependent decay that became considerably more 
rapid than an exponent of — 1 as the starting temperature was 
reduced. 

Here we show that an exact solution to the Fokker-Planck 
equation for instantaneous quenches can explain much of the 
apparently anomalous behavior listed above. The solution 
yields a key insight that the quench dynamics depends cru- 
cially on the vortex-pair distribution function at the initial tem- 
perature Ti. The logarithmic corrections to the vortex decay 
at Ti = Tkt are found to arise from logarithmic corrections 
to the power-law behavior of the initial distributions, which in 
turn comes from the increasing renormalization of the super- 
fluid density near Tkt- 

The starting point is the Fokker-Planck equation for the vor- 
tex pair distribution function 1 14] as used in Ref. 11311 . 
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where T(r, t) is the distribution function for pairs of separa- 
tion r at time t from the instantaneous quench, with r mea- 
sured in units of the vortex core radius ao, t in units of the 
diffusion time a^/2D with D the vortex diffusion constant, 
and r in units <2q. In the limit of the very low vortex den- 
sities that we consider, a term quadratic in T corresponding 
to the recombination of opposite-sign vortices on differing 
vortex pairs 11411 can be neglected in Eq.©. Only the re- 
combination of same-pair vortices at the closest separation 
r = 1 plays a role in the dynamics in this limit. The fac- 
tor K = h 2 a s /m 2 kBT is the dimensionless areal superfluid 
density. The starting value of T at r = 1 and t = is 
To = cxp(— ir 2 Koi/2), using the Villain approximation for 
the vortex core energy. At longer length scales K is renor- 
malized from its initial value Koi by the Kosterlitz scaling 
relation ||7[] 
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At t = the initial condition is that the system is equili- 
brated to a thermal bath at the temperature T,, corresponding 
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to the initial Ki as determined from Eq. ||3}. For quenches 
starting from below Tkt we can make the approximation 
that K i is effectively a constant, since the spatial renormal- 
ization from Eq. ([3]) is rapid, changing from the initial K Qi 
to the renormalized Ki over a length scale that can be less 
than a core radius. The equilibrium distribution then varies 
as T' a r ~ 27TKi , where T' is slightly smaller than Tq due to the 
renormalization. Since the KT transition is a line of critical 
points, all quenches starting at Tj < Tkt are critical quenches 

Immediately after the quench, K in Eq. (O takes on the new 
much larger value Kf corresponding to the low final temper- 
ature Tf. Again, for quenches starting from below Tkt the 
time dependence of Kf is minimal (as well as the space de- 
pendence), recovering after the quench to a superfluid fraction 
of one in less than a diffusion time II 1311 . The thermal bath now 
acts as a delta-function sink to absorb the out-of-equilibrium 
smallest pairs at r = 1 as they annihilate there, which goes 
on until the pair density falls to the equilibrium value at Tf 
(which for a quench to e.g. Tf = 0.1 Tkt can be 20 orders of 
magnitude smaller than the starting value). 

In the limit of constant Ki and Kf, we can solve Eq. (O for 
t > by separation of variables, as detailed in the Supple- 
mental Materials. The solution for the distribution function 
becomes 



r(r,i)=0i-Pi 
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where \Fi is the confluent hypergeometric function of the first 
kind, z = 2 exactly, and 



/ G(l + TvKf - TTK t ) 
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with G the Euler gamma function. This distribution evaluated 
in Fig. la for a quench from T, = 0.9 T KT to T f = 0.1 T KT , 
corresponding to parameter values V = 2.576 x 10~ 4 , Ki — 
0.814, and Kf — 7.479. With no adjustable parameters this is 
seen to give a precise description of the numerical results, and 
the dependence on r/i 1 / 2 exactly satisfies dynamic scaling. 

The time dependence of the vortex pair density (in units o.q) 
is found by integrating the distribution function, 



p{t) = f™T(r,t) 2irrdr 



, q( *(2**,f) 
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The hypergeometric function in this case rises rapidly from 
zero until the time t w 1 where it becomes constant and equal 
to one. Beyond that point the time dependence is then accu- 
rately t~ Zsr - al "l z where 



2nK, 
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is the dynamic exponent first considered by Minnhagen and 
co-workers 01511 . This is just the variation observed earlier 
M 1 311 . but at that time it was not clear that in fact two dynamic 



exponents, z and z sca i e , are involved in the solution. Figure 
lb shows the time decay of the vortex density for the same 
parameters as la, where at long times the decay is i~ 1557 . 

The exact solutions show that Eq. (Q~|) is not correct for 
quenches starting below Tkt- Right at Tkt (and probably at 
higher temperatures), z sca i e takes the value of 2 at long length 
scales, and it is certainly possible that this could be related to 
n = 2 of the superfluid universality class. However, we are 
unaware of any direct relation between these quantities. 

At initial temperatures higher than about 0.95 Tkt the 
above analysis begins to fail because the initial vortex dis- 
tribution is no longer accurately a power law. Right at Tkt 
there is a strong logarithmic correction to Ki such that 2-kKi 
approaches the KT value of 4 (Ki = 2 /if) only at extremely 
large pair separations (it is 4.373 at r = 10 and 4.197 at 
r = 1000). It is still possible to proceed if one notes that 
in Eq. (0 the derivatives are only important near rf»(, since 
at smaller length scales the distribution is nearly flat, while at 
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FIG. 1 . Comparison of the numerical and analytic solutions for the 
vortex distribution function and vortex pair density as functions of 
time, for an instantaneous quench from 0.9 Tkt to 0. 1 Tkt ■ 
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FIG. 2. Comparison of the numerically computed logarithmic slope 
of the vortex decay curve at Tkt in Ref. 1 13] with the exponent 
—Zscaie/z using Ki(r) at the time t = (r/£,o) 2 . Roundoff error 
becomes appreciable at long times in the data since p becomes very 
small. 

larger scales it falls off rapidly. The solutions of Eqs. © and 
© should remain approximately valid, but now 2%Ki has to 
be regarded as the local value at r ~ £, i.e. it becomes time- 
dependent. This in turn means that the vortex density decay 
exponent —z sca i e /z will be time-dependent, approaching —1 
only at very long times. 

We can check this scenario for a quench from Tkt by nu- 
merically evaluating Ki from Eq. |3) as a function of length 
scale and computing the decay exponent — {2%Ki — 2)/2 at 
the corresponding times found from r = £. This can then be 
directly compared to Eq. © by numerically finding the slope 
of the vortex density decay curve for the quench from Tkt 
in Fig.|2]of Ref. II 1 311 as a function of time, shown as the data 
points in Fig. 2 of this paper. Since £o is n °t uniquely deter- 
mined by the above analysis, we have matched the computed 
exponent (solid curve in Fig. 2) to the data at one point, ar- 
bitrarily taken to be t — 100, which yields £o = 11 -5. The 
agreement over more than three orders of magnitude in t be- 
tween the data and the computation shows with little doubt 
that the l/(i In t) behavior seen in [ 1311 was entirely due to the 
log corrections to Ki in the initial distribution. A scaling plot 
of t 2 ^ K *' z T versus r/£ for this case (shown in the Supple- 
mental Materials) is found to give accurate data collapse if £ 
is taken to be ^o^ 1 ^ whereas including a log term in £ does 
not yield data collapse. 

Similar ideas can be applied to quenches from well above 
Tkt, though in that case less is known about the vortex distri- 
bution, and the initial superfluid density is zero. We know at 
Tkt that vortex pairs of extremely large separation (r — > oo) 
distributed at random ("free" vortices) have a pair distribution 
falling off as r~ 4 , and so we expect for temperatures above 
Tkt that this distribution will hold for large but finite r. Fig. [3] 
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FIG. 3. Pair distribution function for a quench from 2.0 to 0.1 Tkt 
(upper plot). The lower plot shows the time dependence of the scale- 
dependent superfluid density. 



shows a numerical evaluation of Eqs. (O and (0) for an instan- 
taneous quench from 2.0 to 0.1 Tkt assuming the initial dis- 
tribution is r~ 4 at all length scales. As the smallest pairs de- 
cay the scale-dependent superfluid fraction quickly recovers at 
the smallest scales, and then with increasing time recovers at 
longer and longer scales as the larger pairs decay. The distri- 
bution function shows a small "hump" at the scales where the 
superfluid fraction is varying most rapidly, which appears very 
similar to the small "humps" seen in the simulations It 1 QTI of 
quenches from 2 Tkt- Integration of these distribution curves 
gives a vortex density decaying accurately as t -1 - for t > 10. 

At short length scales above Tkt there will still be signif- 
icant correlations between the vortices. The pair distribution 
has been found in the simulations of Ref. II 1 QTI at 2 Tkt t° 
fall off as r -3 - 5 for r < 50. Iterating the equations using this 
form in the quench from 2 Tkt gives results entirely similar to 
Fig. 3, but now upon integrating the vortex density falls off as 
t~ o e , just about the behavior seen at short times in the simula- 
tions. We are unable to carry out a more complete calculation, 
however, since we do not know more precisely how the distri- 
bution changes over from r~ 3 - 5 to r~ 4 at large length scales. 
A logarithmic change between these two exponents could well 
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give rise to the \nt/t variation that is observed. Simulations 
on much larger lattice sizes will be needed to study the distri- 
bution function at large length scales. 

The exact solutions presented here for the quench of 2D XY 
superfiuids validate the conclusion in Ref. 11311 that there is no 
"creation" of vortices in a quench, but only monotonic decay 
of the existing thermal vortices. There is no term in the so- 
lution that increases the number of vortices, and which might 
have possibly been missed in the numerical approximations. 
This is contrary to the predictions of the Kibble-Zurek sce- 
nario Biol that vortices will be created in a superfiuid quench 
to low temperatures, with higher densities appearing the faster 
the quench rate. The numerical studies with finite quench 
times 1 13] found just the opposite result, that actually more 
vortices were left over following a slower quench, since in 
that case the system spends more time at higher temperatures 
where the thermal vortex density is higher. The present results 
verify that no excess vortices are created even for an instan- 
taneous quench from above Tkt- The same result has now 
been seen in XY model simulations [10] for quenches start- 
ing from 2 Tkt, where again there was only monotonic de- 
cay of the initial vortex density, and the variation with quench 
rate was entirely similar to the previous numerical results B 1311 . 
The problem in the Kibble-Zurek argument is the restriction 
to measuring the vortex density only at the "freezeout" sam- 
pling time, which increases with the quench time. But since 
the the pairs continually decay, of course this will always re- 
sult in lower vortex densities for a longer quench time and 
hence a later sampling time. But in fact the vortex densities 
can be measured at all times t, as shown in the results above, 
and it then becomes quite clear that the instantaneous super- 
fluid quench has the lowest vortex density at all times of any 
quench rate, since it most rapidly gets to the lowest tempera- 
ture. 

We have also been able to study quenches from low to high 
temperatures Tf < Tkt, at least numerically. The vortex den- 
sity increases initially as i 1 / 2 , indicating proportionality to the 
dynamic length £, and then becomes nearly constant close to 
the equilibrium value at Tf for times greater than a few dif- 
fusion times. The distribution shows the smallest pairs being 
rapidly populated, as expected from the thermal injection at 
r = 1, with the larger pairs only appearing at late times. Fig- 
ures of these results are shown in the Supplemental Materials. 

In summary, we propose the first exact solution for the vor- 
tex dynamics following a critical quench starting from be- 
low Tkt in two-dimensional XY superfiuids. The solution 
highlights the key role of the initial vortex distribution, and 
this allows a consistent explanation of the logarithmic devia- 
tions at Tkt seen in earlier work as being due to logarithmic 
corrections in the initial distribution. It will be important to 
check the validity of the solution by carrying out XY model 
simulations for quenches starting below Tkt, looking for the 
increasingly rapid decay of the vortex density predicted by 
Eq. (fT4l >. Experimentally, these predictions are probably im- 
possible to test in superfiuid helium films where the diffusion 
time is only tens of picoseconds, but may possibly be accessi- 



ble in the 2D superfiuid polariton condensates lfl7ll . which can 
be created and the vortices probed at these rapid time scales. 
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SUPPLEMENTARY MATERIAL 
Exact solution 

In the limit of constant Ki and Kf we can solve Eq. (2) of 
the main text for t > by separation of variables, writing 



T(r,t)=F(w)f(t) 



(8) 



where w — r/t 1 ^ is the dynamic scaling variable. Using the 
chain rule, the Fokker-Planck equation then becomes 



fdt 
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F 
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z t dw 



(9) 

This equation separates if and only if z = 2, and taking the 
separation constant as — a/z yields the two equations 



ldf 
f dt 
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These are straightforward to solve, and the solution for ( fTTT i 
finite at w = (i.e. t — > oo) is 



F(w) = /3 1F1 



, 1 + ttK 
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"27 



(12) 



where the function iF± is the confluent hypergeometric func- 
tion of the first kind, and j3 is a constant. To satisfy the bound- 
ary condition T(r, 0) = T' Q r ~ 27rK ' at t = 0, the large-iu ex- 
pansion of the hypergeometric function requires a = 2irKi , 
and 



P = r 



, / Gjl + irKf-nKi) 
\{2zY K i G(l + TrKf) 



(13) 



where G is the Euler gamma function. The full solution for 
the distribution function then becomes 



irKi, 1 + ttK 



Si 



2z t 2 /* 



f.-2-rrKi/z 

(14) 



which is Eq. (4) of the main text. 



Scaling for quenches from Tkt 

FigureSh shows the scaling plot of t 27rKi / z T versus r/^ for 
an instantaneous quench from Tkt to 0-1 Tkt, using the data 
from Fig. 1 of Ref. [13]. Here £ is taken to have the form £ = 



^t x / z with £o =11.5, and Ki is the value at the scale r w £. 
It is seen that the curves collapse nearly completely, with only 
very slight deviations for r/£ < 1, within the uncertainties of 
the data. Figure |4j3 shows the same plot now using a log form 
for the dynamic scale, £ = £o(t lr^t/to)) 1 , and taking ta = 
0.02. The curves clearly do not satisfy scaling, and the same 
holds true for all values of to < 1- 



CM 




(^ = ^ot 1/Z ) 




1/z. 



r/i;' = ^ (t ln(t/t )) ) 



FTG. 4. Scaling plots of the vortex pair distribution for a quench from 
Tkt- 



Quenches up in temperature 

Figure [5] shows the numerical results for quenches upward 
from 0.1 Tkt to final temperatures of 0.8 and 1.0 Tkt- 
Fig.|5^ shows the growth of the distribution function for the 
quench to 1.0 Tkt, where the smallest pairs are populated 
first from the thermal bath injection at r = 1, and then at 
longer times the larger pairs form to approach the equilibrium 
distribution. Fig.[5j5 shows the resulting time dependence of 
the vortex density, increasing as t 1 ^ 2 , e.g. proportional to the 
dynamic length. The equilibrium thermal density is reached 
within a few diffusion times. 




FIG. 5. Quenches up from 0.1 Tkt to higher temperatures: a) pair 
distribution function for a quench to 1.0 Tkt-, b) vortex density. 



